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Abstract 

Contributions of fermions to the mass of the scalar glueball ++ are calculated at two-loop level 
in the framework of QCD sum rules. It obviously changes the coefficients in the operator product 
expansion (OPE) and shifts the mass of glueball. 
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I. INTRODUCTION 



Quantum Chromo-dynamics (QCD) predicts the existence of glueballs. After a long 
time of experimental and theoretical exploration for glueballs, there is no obvious evidence 
to confirm its existence yet, even though people find several glueball candidates, such as 
/o(1710) and 77(1405) [3] - There are also predictions on the mass of glueballs in various 
theoretical frameworks. Among all the theoretical approaches, the estimation on the mass of 
glueballs by the Lattice QCD and QCD sum rules, seems to be closer to reality. The Lattice 
QCD predicts the mass of the scalar glueball (0++) as 1.5 ~ 1.8GeV @, H, B, H, H, B H, 0, 0- 
With the QCD sum rules, Novikov and Narison evaluated the mass of scalar glueball as 
700 - 900Mev[ll|, Il3fl . whereas Bagan and Steele considered the radiative corrections 
and obtained the mass as 1.7GeV 141] . Later, based on Bagan and Steele's work, Huang 
et al, re-estimated the mass and found a small shift to ~ 1.66GeV|15j. The difference is 
so large that one has reason to doubt if there indeed exists theoretical discrepancy. One 
compelling motivation is that one needs to make a complete calculation which should add 
up the contributions which were neglected in previous calculations. That is the aim of this 
work. 

We have repeated Bagan and Steele's derivations, and noticed that they neglected the con- 
tributions of the fermions to the radiative corrections. Namely, they neglected contributions 
of the loops involving fermion propagators and quark condensates by setting Cf which is a 
coefficient related to quark flavors, to be zero in 14J. They argued, such contributions were 
small compared to others. In this work we include the contributions of the loops containing 
fermion propagators and quark condensates to the correlation function n(g 2 ). With this 
correction, we set a proper platform at sq = 3.8 GeV, where sq represents the threshold for 
the continuum states, and eventually we determine the mass of the glueball as 1721 MeV. 
This value is compatible with that obtained by Bagan and Steele, a bit larger than that 
Huang et al. achieved, but within a tolerable error region, all of them are consistent with 
each other. 



II. SCALAR GLUEBALL QCD SUM RULE 

The correlation function for scalar glueballs is defined as: 

Il( q 2 )= l Jd 4 xe^(0\T{j(x)j(0)}\0), (1) 

where, j(x) = a s G^ u (x)G afJ,u (x) stands for the current of the ++ glueball. By the operator 
product expansion (OPE), the correlation function can be further written as: 

n(g 2 ) = £c»(g 2 )<o|d„|o) 
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[a + 0l log %} (Q 2 ) 2 log %+[b + h log %} (a s G 2 ) 4 



Q 2 ] {g s G 3 ) d 2 
v 2 ^ Q 2 Q 2 



+ Cl log^-]^-^ + -£ (afG 4 ) + ... ; Q 2 = -q 2 > 0, (2) 



where, C n (q 2 ) are the Wilson coefficients, and those operators O n have already well defined 
in ref. 161 ]. A more convenient function form TZk{r, Sq) may be used in later calculations: 

7T J so 

= - H ds s k e- ST ImU(s), (3) 
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where, L the Borel transformation, r is the Borel parameter, and Sq represents the threshold 



n k (r,s ) = -L[(Q 2 ) k Il(-Q 2 )]-- / dss k e- ST ImU(s) 

T 7T Jso 



for the continuum states, substituting Eq ([2]) into Eq Q, then we have 14| (for k > — 1): 



ft_i(T,a ) = -^[l-pi(s r)] + 2^{ 7i; + J E 1 ( So r) + log So r + e-^ 

- 1 - [1 - Pl ( So r)] log ^} + n(0) - {6 - 6i[7s + logrz/ 2 

+ E x (s r)} } (a s G 2 ) - {c + Cl [l- lE - log rv 2 - E x (s r) (4) 

S T 2 

K k (r,s ) = (-jt^TU^so) (5) 

OT 

h 

where, p fc (x) = e" x ^ ^ J /j\, E^x) = / dye~ y /y, (6) 

3=0 ^ 

and, 7b = Euler's constant ~ 0.5772. (7) 



With the function TZk(r, sq), the mass of the scalar glueball is: 

9 n t> \ ^fe+l(T,S ) 



MqJt, s ) = -—(\o g n 



Although the mass of glueball is determined by a sum of 7?.^ with different integer k values, 
only TZi is mostly important and kept in the final result. The reason is that only IZx is 
reliable for determination of the ++ glueball mass[l5|. According to this comment above, 
for the 1Z with k > 0, only a, (i = 0, 1), b\ and C\ in Q can affect the final mass of the 
glueball, since other coefficients will vanish through the derivation in EqfJHJ). 

III. QUARK CONTRIBUTION IN CORRELATION FUNCTION I% 2 ) 

With the fermion contributions containing the fermion propagators and the quark con- 
densations, the correlation function II (q 2 ) can be divided as: 

n(g 2 ) = nV) + nV) (9) 
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(a) (b) (c) 

FIG. 1: (a) the pertubative part; (b) the quark condensation part; (c) the gluon condensation part 



where, the II 1 (q 2 ) is the correlation function without fermion contributions, and was given 
by Bagan and Steele 0- The Uf(q 2 ) stands for the fermion part. With OPE, we have 

ipV) = c{{d ) + c((d 3 ) + c{(6 A ) + (io) 

where, O n is defined above. Cq is the Wilson coefficient for the unit operator Oo in Fig 1-a, 
and in MS scheme, we have: 

V2 



CM) = (^) 3 (Q 2 ) 2 (51og%-21og 2 %). (11) 

7X V* V A 

Cl is the Wilson coefficient for operators with quark condensate. From Fig(l-b), we have: 



-^[-m^mjr) - E 1 (s r) + logm 2 (l - p (0))) + -^( P o(rs ) 
Po«r)) + ^4%(p 2 Kr) - l)](qq) + ^[^1(- El (m 2 q r) + E^rso) 
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2 40 
- logm 2 (l - p (0))) + — (ftKr) - 1) + 5-5-3(1 - p 2 (m 2 r))] 
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x (g s qqG) + a 2 [— (-(^i(mjr) - #l(ts )) - logm 2 (l - p (0))) 

+ ^^( 1 -^K^)) + ^p(^K^)-i)](^> ) (12) 

where, pk(x) and ^(x) are defined in Eq and g stands for the u, d and s quarks. Later, 
we will show that, the contributions from the parts with quark condensates are less than 
1%, so in general, can be safely ignored. 

C((q 2 ) is the Wilson coefficient for the operators with gluon condensates. Generally, there 
are two ways to calculate the Wilson coefficients (l7| . one is the plane wave method, by which 
Il 1 (g 2 ) is given (See Eq Q). The another is the fixed-point gauge technique, by which we 
have: 



n, 2 13 — n 2 

?/^//Y\ - ^(_i£ + 2 log-^ 

7T 3 V 1 



Cl(q 2 )(0 4 ) = -JL(- + 21og-^)(a s G 2 >. (13) 
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FIG. 2: (a) the mass without C( at s = 3.8GeV 2 ; (b) the mass with C[ at s = 3.8 GeV 2 



Substituting all the corrections ( (ITT]) . (IT2"1) . and (TT3"j) ) back into Eq ([2]), we have the correlation 
function IT(g 2 ) and the various functions TZk with corrections from fermions. The result is 
shown in fig(2): We choose the reasonable platform at r G {0.4, 0.6} in the region 3.6 GeV 2 < 
So < 4.2 GeV 2 . Within the platform, the ratio of the contribution of the unit operator term 
Co(g 2 )Oo, which stands for the pertubative part, to the mass determined by IZk is more than 
90%, it enables the OPE expansion to converge sufficiently fast. Besides, in the region, the 
ratio of the contribution of the continuum part to the mass is less than 30% — 40%. It implies 
that this value of sq is appropriate for the quark-hardron duality. Within this platform we 
have obtained the ++ glueball mass as: 1.721 ± 0.065 GeV 2 , where the error is caused by 
the variable of the Sq in the region, meanwhile, the error caused by the variable of the Borel 
parameter r is very tiny so that it can be ignored. Fig (2-a) and Fig (2-b) show that, the 
quark condensate contributes little to the mass of the glueball, since the C{ given in Eq (fT2l 
turns to zero at mq — > 0. That is why we could directly determine the mass of glueball by 
neglecting C3 as long as the mass of light quarks is small. 



IV. CONCLUSION AND DISCUSSION 



In this paper, we analyze the contribution of the diagrams involving internal fermion lines 
and quark condensates. Following the traditional way, we determinate the mass of the ++ 
scalar glueball as 1.728 ± 0.132 GeV 2 . Comparing the result of Huang (m ~ 1.66GeV 2 ), a 
little shift of the mass is resulted in by taking the fermion condensation into accout. 
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